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The production of K® mesons in pp collisions at a centre-of-mass energy of 0.9 TeV is 
studied with the LHCb detector at the Large Hadron Collider. The luminosity of the anal- 
ysed sample is determined using a novel technique, involving measurements of the beam 
currents, sizes and positions, and is found to be 6.8 ±1.0 /ib~ . The differential prompt 
-Kg production cross-section is measured as a function of the Kg transverse momentum 
and rapidity in the region < p^ < 1.6 GeV/c and 2.5 < y < 4.0. The data are found to 
be in reasonable agreement with previous measurements and generator expectations. 
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1 Introduction 



Strangeness production studies provide sensitive tests of soft hadronic interactions, as the 
mass of the strange quark is of the order of Aqcd- Strange-hadron production is sup- 
pressed, as a consequence, but still occurs in the non-perturbative regime. The hadronic 
production of K$ mesons has been studied by several experiments at a range of different 
centre-of-mass energies, both in pp and pp collisions (see for example (lH2|). The most 
recent measurements of K$ production at the Tevatron have shown deviations with re- 
spect to the expectations of hadronization models 0. Strangeness production is also a 
topic of great interest in heavy ion physics, and measurements of this process in pp and 
pp collisions serve as reference point [7J. 

In this paper measurements of prompt Kg production are presented using data col- 
lected with the LHCb detector in pp collisions at y/s = 0.9 TeV, during the 2009 pilot 
run of the Large Hadron Collider (LHC). A K® is defined to be prompt if it is directly 
produced in the pp collision, or if it appears in the decay chain of a non-weakly-decaying 
resonance (such as K*) directly produced in the pp collision. The measurements are made 
in the rapidity interval 2.5 < y < 4.0 and down to below 0.2 GeV/c transverse momentum 
with respect to the beam line. This is a region not explored at this energy by any previ- 
ous experiment, and is complementary to the coverage of other LHC experiments. The 
determination of the prompt K§ production cross-section is normalized using an absolute 
measurement of the luminosity that relies on knowledge of the beam profiles. 

The paper is organized as follows. Section [2] gives a brief description of the LHCb 
detector and the configuration used to record data in December 2009 during the LHC 
pilot run. Section |3] gives an overview of the analysis strategy, the details of which are 
presented in the three following sections. Section H] is dedicated to an explanation of the 
luminosity measurement, Section [5] presents the Kg candidate selection and Section 
the determination of the Kg trigger and reconstruction efficiencies. The final results 
are discussed in Section [7] and compared with model expectations, before concluding in 
Section |HJ 



2 LHCb detector and 2009 data sample 

The LHCb detector is a single-arm magnetic dipole spectrometer with a polar angular 
coverage with respect to the beam line of approximately 15 to 300 mrad in the horizontal 
bending plane, and 15 to 250 mrad in the vertical non-bending plane. The detector is 
described in detail elsewhere |. All subdetectors were fully operational and in a stable 
condition for the data that are analysed. For the measurements presented in this paper 
the tracking detectors and trigger strategy are of particular importance. 

A right-handed coordinate system is defined with its origin at the nominal pp inter- 
action point, the z axis along the beam line and pointing towards the magnet, and the y 
axis pointing upwards. Beam-1 (beam-2) travels in the direction of positive (negative) z. 

The LHCb tracking system consists of the Vertex Locator (VELO) surrounding the 
pp interaction region, a tracking station (TT) upstream of the dipole magnet, and three 
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tracking stations (T1-T3) downstream of the magnet. Particles traversing from the in- 
teraction region to the downstream tracking stations experience a bending-field integral 
of 3.7 Tm on average. 

The VELO consists of silicon microstrip modules, providing a measure of the radial 
and azimuthal coordinates, r and 0, distributed in 23 stations arranged along the beam 
direction. The first two stations at the most upstream z positions are instrumented to 
provide information on the number of visible interactions in the detector at the first level 
of the trigger ('pile-up detector'). The VELO is constructed in two halves (left and right), 
movable in the x and y directions so that it can be centred on the beam. During stable 
beam conditions the two halves are located at their nominal closed position, with active 
silicon at 8 mm from the beams, providing full azimuthal coverage. During injection and 
beam adjustments the two halves are moved apart horizontally to a retracted position 
away from the beams. 

The TT station also uses silicon microstrip technology. The downstream tracking 
stations T1-T3 have silicon microstrips in the region close to the beam pipe (Inner Tracker, 
IT), whereas straw tubes are employed in the outer region (Outer Tracker, OT). 

During the 2009 run, low intensity beams collided in LHCb at the LHC injection 
energy, corresponding to a total energy of 0.9 TeV. Due to the dipole magnetic field the 
beams have a crossing angle that results in the pp centre-of-mass frame moving with 
velocity 0.0021c in the — x direction. Both the beam sizes and crossing angle were larger 
than those designed for high-energy collisions. In order not to risk the safety of the 
VELO, the 2009 data were recorded with the two VELO halves positioned 15 mm away 
from their nominal data-taking position (VELO partially open), resulting in a reduced 
azimuthal coverage. For this run, the magnetic dipole field was pointing downwards. 

The bulk of the data presented here were collected in a series of LHC fills with the 
following two sets of beam conditions. The first configuration contained four bunches per 
beam, spaced by more than 8 /is, with two colliding and two non-colliding bunches, and a 
total peak beam intensity of about 1.8 x 10 10 protons per bunch. The second configuration 
contained 16 bunches per beam, spaced by more than 2 /xs, with eight colliding and eight 
non-colliding bunches, and a total peak beam intensity of about 1.3 x 10 10 protons per 
bunch. The nominal LHC injection optical function at the interaction point was used 
{P* = 10 m). 

A trigger strategy was deployed to provide high efficiency for pp inelastic interactions 
and for beam collisions with the residual gas in the vacuum chamber. The latter class 
of events is a necessary ingredient for the luminosity analysis. Events were collected for 
three bunch-crossing types: two colliding bunches {bb), beam-1 bunch with no beam-2 
bunch (bl), and beam-2 bunch with no beam-1 bunch (62). The first two categories of 
crossings, which produce particles in the forward (+z) direction, were triggered using 
calorimeter information: a 2x2 cluster with more than 240 MeV of transverse energy in 
the Hadron Calorimeter (HCAL) and at least three hits in the 6016 cells of the Scintillator 
Pad Detector (SPD) at the entrance to the calorimeter were required. Events containing 
a track in the muon system with transverse momentum above 480 MeV/c were also 
triggered. Crossings of the type 62, which produce particles in the backward direction 
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only, were triggered by demanding a hit multiplicity of more than seven in the pile-up 
detector. 

The visible collision rate for a single bunch pair was about 10 Hz and the acquired bl 
(b2) rate for a single bunch was approximately 0.015 Hz (0.002 Hz), in agreement with the 
measured residual pressure and VELO acceptance. A sample of 424 193 events triggered 
in bb crossings is used in the K$ analysis. 

3 Analysis strategy 

All K$ candidates are reconstructed in the 7r + 7r~ decay mode, using only events triggered 
by the calorimeter. Contributions from secondary interactions in the detector material 
or from the decay of long-lived particles are suppressed by requiring the K$ candidates 
to point back to the pp-collision point. No attempt is made to separate the contributions 
from K§ mesons produced in diffractive and non-diffractive processes. 

Due to the long K® lifetime and partially open VELO position, only a small fraction of 
the -Kg daughter tracks traversing the spectrometer leave a signal in the VELO. Therefore, 
two paths are followed for the K® reconstruction and selection: 

a) Downstream-track selection: 

Tracks reconstructed only with hits in the TT and T1-T3 stations (called down- 
stream tracks) are combined, without using the VELO. The origin of the K§ is taken 
as the point on the z axis that is closest to the reconstructed flight vector of the 
J<"g candidate. This point is taken as an estimate of the primary vertex (PV), and 
is referred to as the 'pseudo-PV. 

b) Long-track selection: 

K§ candidates are formed with tracks leaving hits in the VELO and in the T stations 
(called long tracks). If available, measurements in the TT are added to the tracks. 
The PV is reconstructed from tracks seen in the detector, using VELO information 
whenever available. 

The analysis is performed in bins of K§ phase space. The kinematic variables used are 
the K$ transverse momentum pr = \/pl + Vy an d the rapidity y = \ \n((E+p z ) / (E—p z )), 
where (E,p) is the K® four-momentum in the pp centre-of-mass system. For a given bin 
i in px and y, the prompt K$ production cross-section is calculated as 

yyobs 

*i = , (!) 

trig/sel ad T ' ^ ' 

e i e i ^int 

where N° hs is the number of observed K$ — > 7r + 7r" signal decays with reconstructed 
and y in bin i, ef 51 the reconstruction and selection efficiency, e* ng / sel the trigger efficiency 
on selected events, and L int the integrated luminosity. The number of signal events iV° bs 
is obtained from the mass distributions of the K§ candidates. 
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The reconstruction and selection efficiency is estimated from a fully-simulated Monte 
Carlo (MC) sample of single pp collisions as 



N scl 

e sel - 1 o) 

i jyP rom Pt ' ' 

where Nf el is the number of K$ — > tt + tt^ signal decays selected in the untriggered MC 
sample with reconstructed pi and y in bin i (extracted using the same procedure as in the 
data), and where ]\rP roinpt j s the number of generated prompt Kg mesons with generated 
Pt and y in bin i. This efficiency includes the geometrical acceptance, as well as the 
reconstruction and selection efficiencies. It also incorporates all corrections related to the 
following effects: secondary interactions of K§ in the material, K® — > 7T + 7T~ branching 
fraction, decay in flight and secondary interaction of the decay products, non-prompt K<$ 
production and finite resolution of the px and y observables. 

The trigger efficiency is estimated using the same MC events. However, since the 
efficiency depends on the global event properties, the MC events are weighted to reproduce 
the observed track multiplicity in the selected signal events. Then 

trig/scl _ *i m\ 

€ i - vscl ^ ' 

i 

is computed, where y. tng / scl and y^ sel are the weighted MC signal yields extracted after 
and before the trigger cuts are applied. 

The integrated luminosity Lj nt is determined using a novel 'beam imaging' method jij, 
taking advantage of proton collisions with the residual gas in the interaction region and 
of the excellent vertexing capability of the VELO. The beam profiles and positions are 
reconstructed using tracks produced in beam-gas and beam-beam collisions. Combining 
this information with bunch current measurements from the LHC machine yields a direct 
measurement of the integrated luminosity. 



4 Luminosity determination 

In the relativistic approximation, the average instantaneous luminosity produced by one 



pair of colliding bunches can be expressed as [10 



L = 2 c ni n 2 f cos* »/„(«,»,*, « W«, V,*,t)i* d» d; it , (4) 

where rij are the number of protons in bunch % (i = 1,2), / = 11.245 kHz is the LHC 
revolution frequency, 9 is the half crossing angle of the beams, and pi(x,y,z,t) is the 
density of bunch % normalized as J pi(x,y,z,t)dxdydz = 1 at all times t. The overlap 
integral in Eq. (j4j) is taken over the duration of one bunch crossing. Tracks measured in 
the VELO allow vertices from beam-gas and beam-beam collisions to be reconstructed 
for each pair of bunches. From the distributions of these vertices, and assuming the gas 
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Figure 1: Distributions in the horizontal (top) and vertical (bottom) planes of the reconstructed 
vertices in 61 (blue filled circles and solid fit line) and b2 (red open circles and dashed fit line) 
crossings in one fill. 



density to be uniform in any plane transverse to the beams, the positions, angles and 
sizes of the bunches are measured, and their overlap integral is computed. The numbers 
of protons per bunch are determined with the LHC machine instrumentation, enabling 
an absolute normalization of the luminosity. The total luminosity is then obtained as the 
sum of the estimates for each pair of colliding bunches in the machine. 

The beam crossing angle is limited to the horizontal plane. No correlation between 
the transverse coordinates is observed at the level of precision needed for this analysis, 
thus the x and y projections can be factorized. The bunch shapes are well described by 
Gaussian distributions in all three dimensions, characterized in the x — y plane at the time 
of crossing by their width and their mean position ^ (j = x,y), and by their average 
longitudinal width a z , assumed to be equal for both beams. With these approximations 
and for small crossing angle, Eq. (jlj) can be rewritten as 

L = ^2 f tt 1 exv (_}i ^-^y 



The observables cx^ and fly are extracted from the transverse distributions of the beam- 
gas vertices reconstructed in the bb crossings of the colliding bunch pair with a z coordinate 
satisfying -1000 < z < -200 mm (200 < z < 1000 mm) for % = 1 (i = 2). These 
transverse distributions are obtained by projecting the reconstructed vertex positions 
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Table 1: Parameters describing the vertex resolution functions defined in Eqs. (|7|) and ([8]). The 
quoted errors include statistical and systematic uncertainties. The parameters fj and rj were 
fixed in the fits, and hence have no uncertainties. 







r j 


s track 


Si 




mij [m *] 


b 2j 


m 2j 


m 1 




X 


0.9 


0.32 


177 ±7 


5.9 ± 1.1 


1.18 ±0.07 


-0.86 ±0.30 


0.83 ±0.14 


+0.77 ±0.24 


y 


0.9 


0.36 


164 ±6 


3.7 ±1.1 


1.24 ±0.08 


-0.57 ±0.16 


0.85 ±0.14 


+0.77 ±0.24 



onto a plane perpendicular to the corresponding beam direction. As illustrated in Fig. [TJ 
the beam directions, and hence also the half crossing angle 9, are obtained from straight- 
line fits through the measured positions of vertices reconstructed in bl and 62 crossings 
of other non-colliding bunches. The observed half crossing angle of 6 = 2. 1 ± 0.1 mrad in 
the horizontal plane is in agreement with the expected value. 

In addition, the distribution of pp-collision vertices, produced by the colliding bunch 
pair and identified by requiring —150 < z < 150 mm, can be used to measure the 
parameters of the luminous region. Its position \ij and transverse width <jj, 

& = ^ll' T ^T 13 and a j = ' U = x,y) (6) 

can be used to constrain the bunch observables. Owing to the higher statistics of pp 
collisions compared to beam-gas interactions, the constraints of Eq. ([6]) provide the most 
significant input to the overlap integral. 

The longitudinal bunch size a z is extracted from the longitudinal distribution of the 
pp-collision vertices. Because a z is approximately 200 times larger than cr^, the crossing 
angle reduces the luminosity by a non-negligible factor equal to the first square root term 
in Eq. (J3]). For the fill used to determine the absolute luminosity, this factor is estimated 
to be 1.087 + 0.012. 

The vertex resolutions need to be measured since they are of the same order as the 
bunch sizes. This is achieved by comparing, on an event-by-event basis, the reconstructed 
vertex positions obtained from two independent sets of tracks. In each event, the sample of 
available tracks is randomly split into two sets of equal multiplicity, and the event is kept 
only if exactly one vertex is reconstructed for each set. In this case the two vertices are 
assumed to originate from the same interaction. The vertex resolution for each coordinate 
is obtained as the width of the distribution of the difference in position between the two 
reconstructed vertices divided by \/2- A systematic study of the vertex resolutions in 
both x and y is then performed as a function of the number of tracks N contributing to 
the vertex, of the crossing type, and of the z coordinate of the vertex. The resolution 
functions are found to be well parametrized by a double Gaussian function 

Rj(N, z) = / ■ G(sj(N, z)) + (1 - fj) G( Sj (N, z)/ rj ) , (j =x,y), (7) 

where fj is the fraction of events in the first Gaussian function, rj is the ratio of the 
widths of the two Gaussian functions, and G(sj(N, z)) is a Gaussian function centred at 
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Figure 2: Measured profiles of beam-1, beam-2 and luminous region (from left to right) in the 
horizontal direction x, in bb crossings of one pair of colliding bunches in one fill. The solid curve 
is a fit to the observed distributions, the dotted curve represents the vertex resolution, and the 
dashed curve shows the underlying distributions after deconvolution of the vertex resolution. 



zero with width 

8 f(N, z) = AT-o.s+^/tv 2 s track for beam-beam .. _ . . . 

s)(N,z) = (hj + m ijZ ) sf(N,z) for beam-gas (z = 1,2) ' v- x >y)- 

The parameters s*. rack describe the per-track resolutions, Sj specify the dependence on the 
number of tracks, while bij and m^- model the linear z dependence for beam-gas vertices. 
The validity of this parametrization has been verified with MC simulation studies. The 
systematic uncertainties on the parameters are estimated from the level of agreement in 
that check. The final set of resolution parameters is given in Table [TJ The resolution 
is found to be better in y than in x, which is expected from the partial VELO opening 
described in Section [2j 

For both transverse coordinates, each sample of vertices (defined for every colliding 
bunch pair in each fill) is fitted with convolutions of the Gaussian beam shapes with the 
resolution function of Eq. (J2J). This fit is performed with all three types of interactions. 
With the constraints of Eq. dHJ), this yields directly the position and Gaussian width 
o\j of the underlying distributions. Some example distributions are shown in Fig. [2j The 
systematic errors on the results are estimated by varying the resolution parameters within 
their total uncertainties. 

The remaining ingredients needed for the direct luminosity measurement are the 
bunch intensities. The LHC is equipped with two systems of beam current transformers 
(BCT) (llj . A DC-BCT system provides an ungated measurement of the total beam 
current, while a fast-BCT system is gated to measure the current induced on a bunch- by- 
bunch basis. The individual bunch intensities are obtained from these fast-BCT readings, 
but constraining their sum to the DC-BCT measurements. At the low intensities of the 
2009 pilot run, the offset in the DC-BCT digitization is non-negligible and is corrected 
by averaging the readings in the periods without circulating beam just before and after a 
fill. 
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Table 2: Kg — > 7t + -k selection requirements. 


Variable 


Requirement 


Downstream-track selection 


Each 7T-track momentum 

Each 7r-track transverse momentum 

Each track fit x 2 /ndf 

Distance of closest approach of each 7r-track to the z axis 
Kg decay vertex fit x 2 /ndf 
z of Kg decay vertex 
\z\ of pseudo-PV 

COS ^pointing 

Kg proper time (cr) 


> 2 GeV/c 

> 0.05 GeV/c 

< 25 

> 3 mm 

< 25 

< 2200 mm 

< 150 mm 

> 0.99995 

> 5 mm 


Long-track selection 


z of associated PV 
Each track fit x 2 / n df 
Kg decay vertex % 2 /ndf 
z(K$) - *(FV) 

Variable v related to impact parameters 


< 200 mm 

< 25 

< 100 

> mm 

> 2 



The method described above was used to measure the luminosity in four different 
machine fills. Two of those fills were relatively short and the third was taken before 
optimization of the beam alignment. The remaining fill, taken under optimal conditions 
and representing approximately 25% of the sample used for the Kg production study, is 
chosen to determine the absolute normalization of the luminosity for the data set used 
for the Kg analysis. The other three fills yield less precise but consistent results. The 
integrated luminosity for the data set used for the Kg selection, L int = 6.8 ± 1.0 /zb _1 , is 
obtained by scaling with the number of pp interaction vertices measured with the VELO. 
The relative uncertainty on this result comprises contributions from the measurements of 
the beam intensities (12%), widths (5%), relative positions (3%) and crossing angle (1%). 
This is the most precise determination of the luminosity for the 2009 LHC pilot run. The 
limiting uncertainty on the beam intensity is expected to improve in the future. 

5 K$ selection and signal extraction 

In the downstream-track selection, a Kg candidate is formed from any combination of two 
oppositely-charged downstream tracks, assumed to be pions, satisfying the requirements 
listed in the top part of Table [2j The pseudo-PV was defined in Section [31 and pointing 
is the angle between the Kg momentum vector and the direction joining the pseudo-PV 
and the Kg decay vertex. 
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Figure 3: Mass distributions of all selected candidates, in the downstream-track (left) and 
long-track (right) selections. The points are the beam-gas subtracted data and the curves are 
the result of the fits described in the text. 



In the long-track selection, primary vertices are reconstructed from at least three 
tracks. Each K$ candidate formed from long tracks is associated with the PV that 
minimizes its impact parameter and the requirements listed in the bottom part of Table [2] 
are applied. The variable v is similar to a Fisher discriminant formed with the logarithms 
of the impact parameters; it is defined as v = In [(J + J_)/(/ Ii)}. Here I + , J_ and 1$ are 
the impact parameters of each of the two tracks and of the K$ candidate with respect to 
their closest PV, respectively, and the value of I\ is fixed to 1 mm. 

Mass distributions are obtained for both bb crossings and 61 crossings. In order to keep 
only the contribution arising from pp collisions, the 61 mass distribution is subtracted, 
after proper normalization, from the bb mass distribution. The normalization factor is 
0.908 ± 0.015, averaged over the entire sample used for this analysis. It is obtained from 
the ratio of the number of interaction vertices reconstructed in bb and bl crossings in the 
region z < —200 mm where no pp collision can take place. This beam-gas subtraction 
removes about 1.2% of the K$ signal. 

The beam-gas subtracted mass distributions are shown in Fig. [3] for all selected K$ 
candidates. A \ 2 fit is made, describing the background with a linear function and the 
signal with the sum of two Gaussian functions of common mean value, with all param- 
eters left free. It gives a total K$ signal yield of 4801 ± 84 (1140 ± 35), a mean mass 
value of 497.12 ± 0.14 MeV/c 2 (497.43 ± 0.14 MeV/c 2 ), and an average resolution of 
9.2 MeV/c 2 (5.5 MeV/c 2 ) for the downstream-track (long-track) selection. Quoted un- 
certainties are statistical only. The mass values are close to the known Kg mass value of 
497.61 ± 0.02 MeV/c 2 Q, reflecting the current status of the mass-scale calibration. In 
the long-track selection, the statistics are lower than in the downstream-track selection, 
but the background level is lower and the mass resolution is significantly better. 

The beam-gas subtraction and signal yield extraction are repeated for each bin in 
Pt and y, leading to the results shown in Table [3j The systematic uncertainties on the 
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Table 3: Number of observed beam-gas subtracted Kg — > ir + ir~ signal decays, as extracted 
in the downstream- and long-track selections for each bin of transverse momentum and 
rapidity y. The first quoted uncertainty is statistical and the second systematic. The latter is 
uncorrelated across bins. A dash indicates that the statistics were insufficient to determine a 
result in that bin. 



p T [GeV/c] 


2.5 < 


y<3.0 


3.0 < y < 3.5 


3.5 < 


y < 4.0 


Downstream-track selection 


0.0- 


0.2 






73 ± 10 ± 2 


40 ± 


8 ± 


1 


0.2- 


0.4 






278 ± 21 ± 6 


288 ± 21 ± 10 


0.4- 


0.6 


147 ± 15 ± 4 


428 ± 24 ± 7 


388 ± 21 ± 10 


0.6- 


0.8 


202 ± 16 ± 1 


379 ± 22 ± 8 


332 ± 21 ± 


8 


0.8- 


1.0 


176 ± 15 ± 1 


213 ± 16 ± 6 


217 ± 17 ± 


1 


1.0- 


1.2 


113 ± 11 ± 1 


173 ± 14 ± 1 


111 ± 12 ± 


4 


1.2 - 


1.4 


94 ± 11 ± 2 


90 ± 10 ± 


32 ± 


8 ± 





1.4- 


1.6 


56 ± 


8 ± 2 


64 ± 8 ± 3 


20 ± 


5 ± 


1 


Long-track selection 


0.0- 


0.2 


17 ± 


5 ± 2 


34 ± 7 ± 3 








0.2- 


0.4 


31 ± 


6 ± 2 


75 ± 9 ± 4 








0.4- 


0.6 


63 ± 


8 ± 6 


121 ± 12 ± 3 


41 ± 


7 ± 


1 


0.6- 


0.8 


64 ± 


8 ± 2 


134 ± 12 ± 3 


65 ± 


9 ± 


5 


0.8- 


1.0 


50 ± 


7 ± 2 


91 ± 10 ± 2 


53 ± 


8 ± 


4 


1.0- 


1.2 


30 ± 


6 ± 1 


40 ± 7 ± 5 


35 ± 


7 ± 


2 


1.2 - 


1.4 


16 ± 


4 ± 


33 ± 6 ± 5 


27 ± 


5 ± 


6 


1.4- 


1.6 


8 ± 


3 ± 


19 ± 5 ± 3 


14 ± 


4 ± 


2 



extraction of these yields are obtained by comparing the yields from single and double 
Gaussian signal fits and from side-band subtraction to the expected yield in a Monte 
Carlo sample of comparable statistics to the data set. Additionally the fitted and side- 
band subtracted yields are compared, and an alternate (exponential) background model 
is used in the mass fits. The largest observed deviation in any of these studies is taken as 
systematic uncertainty. For the long-track selection, where the yields are lower, the central 
value is obtained from the side-band subtraction method assuming a linear background. 



6 Efficiency estimation 

A sample of fully simulated events is used to estimate the reconstruction and selec- 
tion efficiency e^ el in each and y bin. Single pp collisions are generated with the 
PYTHIA 6.4 program [l3[ and the generated particles are tracked through the detector 



with the GEANT 4 package [14J, taking into account the details of the geometry and 



material composition of the detector. The simulation of the detector response is tuned 
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to reproduce test beam results [8|. In terms of dead and noisy channels, the simulation 
reflects the detector status of the data set used in this analysis. 

Residual misalignments of the tracking stations and edge-effects of cell efficiencies in 
the Outer Tracker are not perfectly described in the MC sample, resulting in an over- 
estimation of the tracking efficiency. To incorporate these effects, we compare for each 
detector unit the hit content of the tracks in the data and MC samples and randomly 
remove hits in the simulation until we achieve agreement in all subdetector components 
and phase-space regions. The MC sample modified in this way is the nominal MC sample, 
used throughout the analysis. 

To assign systematic uncertainties on the efficiencies obtained in this MC sample the 
single track-finding efficiencies were measured. The VELO efficiency is obtained by using 
reconstructed tracks in the TT and in the T1-T3 stations and checking for an associated 
track segment in the VELO. Similarly the TT and T1-T3 station efficiencies are tested by 
reconstructing tracks using VELO and HCAL information. For downstream tracks with a 
Pr larger than 0.2 GeV/c agreement between the track-finding efficiencies in data and in 
the Monte Carlo sample is observed within the statistical uncertainties of approximately 
3%. Below 0.2 GeV/c, the ratio of efficiencies in data and MC is found to be 0.85±0.12. As 
a conservative approach 3% (15%) uncertainties for the reconstruction efficiency of tracks 
with a pr larger (smaller) than 0.2 GeV/c are assigned. Propagating these uncertainties 
to the Kg reconstruction efficiency results in correlated systematic uncertainties of up to 
17% for the lowest K§ bins of the downstream-track selection. 

The systematic uncertainty on the K§ selection efficiency is obtained by comparing, 
in data and MC, the selection efficiency relative to a preselection. This preselection is 
close to 90% efficient for downstream-track selected signal events in MC. 

If the reconstruction and selection efficiency varies strongly within a given bin of phase 
space, the average value estimated with MC will depend on the assumed production 
spectrum within the bin. The extraction of the efficiency-corrected yield in each bin is 
therefore repeated using efficiencies in four sub-bins rather than an average efficiency, and 
the difference with respect to the nominal result is taken as an uncorrelated systematic 
uncertainty. The size of this effect varies between and 20%. The largest uncertainties 
are obtained in bins at the limit of the acceptance. 

The fraction of non-prompt K§ signal in the selected MC sample is found to be 0.6%. 
By definition, this is corrected for in the efficiencies defined in Eq. (T2]). Because the 
correction is so small, even doubling this contribution would have no significant impact 
on the final result. Similarly, the systematic uncertainty due to material interactions, 
assuming a conservative ±10% variation of the known detector material, is found to be 
negligible. 

The trigger efficiency e* ng//&el for selected signal events depends on the track multiplicity. 
As outlined in Section 3, e* ng//sel is obtained after weighting the previously-defined nominal 
MC sample in order to reproduce, in selected signal events, the track multiplicity observed 
in the data (see Fig. 2] (left)). This re- weighting is only applied for the determination of 
the trigger efficiency, as the reconstruction and selection efficiency has been shown not to 
depend on the track multiplicity. The trigger efficiency is found to be greater than 95% 
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Figure 4: Left: Downstream track multiplicity for events containing a signal in data (black 
filled circles) and MC (red open squares), normalized to unit area. Right: Trigger efficiency for 
events containing a signal K§ decay in the downstream-track selection, as a function of the 
Pt, estimated both in data (black filled circles) and MC (red open squares), using Eq. ([9]). 



in every phase-space bin. As a cross check, it is also extracted directly from data, using 
a method that exploits the fact that signal events can be triggered by the K$ daughters 
(trigger on signal, TOS) or by the rest of the event (trigger independent of signal, TIS), 
with a very large overlap between the two cases. Assuming that the two ways to trigger 
are independent, A"tis&tos = £tis ctos A sd = A r Tis N TOS /N Beh where, in a given region 
of phase space, Axis and Atos are the number of TIS and TOS events, Atis&tos is the 
number of events which are simultaneously both TIS and TOS, and A se i is the number of 
selected signal events. Hence 

^trig/sel _ A T i S |TOS _ A T IS|TOS A T i S &TOS /^s 

A S ci A TIS A TO s 

where Atis|tos is the number of events which are triggered either as TIS or TOS. Due 
to the limited data statistics, a significant comparison between data and MC can only be 
done in bands of or y, rather than in 2-dimensional bins. Good agreement is found, 
as illustrated in Fig. 0] (right), and the observed differences are translated into a global 
correlated systematic uncertainty of 2%. 

The dependence on the modeling of diffractive processes is studied per bin of phase 
space by changing the fraction of diffractive events in the PYTHIA 6.4 sample by 
50% of its value, and by replacing these events with diffractive events generated with 
PYTHIA 8.1 00. The evaluation of the MC efficiencies is repeated for different 



PYTHIA 6.4 parameter values [17j], leading to no significant change. 

There are two important differences in the analysis of the K$ candidates from the 
long-track selection relative to the downstream-track selection. Firstly, a reconstructed 



2 We consider single- and double-diffractive process types 92-94 in PYTHIA 6.421, which includes only 
soft diffraction, and 103-105 in PYTHIA 8.130 (soft and hard diffraction). 
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Table 4: Total efficiencies (in %) in bins of transverse momentum and rapidity y for the two 
selections. The first uncertainty is uncorrelated, including the statistical uncertainty from MC, 
and the second is at least partially correlated across bins. 



p T [GeV/c] 


2.5 < y < 3.0 


3.0 < y < 3.5 


3.5 < y < 4.0 


Downstream-track selection 


0.0 


-0.2 




3.4 ± 0.5 ± 0.5 


3.0 ± 0.6 ± 0.3 


0.2 


-0.4 




7.3 ± 0.2 ± 0.8 


7.4 ± 0.2 ± 1.0 


0.4 


-0.6 


3.5 ± 0.4 ± 0.4 


11.8 ± 0.2 ± 0.9 


12.0 ± 0.2 ± 0.9 


0.6 


-0.8 


7.4 ± 0.3 ± 0.5 


15.0 ± 0.2 ± 1.2 


15.1 ± 0.2 ± 1.2 


0.8 


- 1.0 


11.1 ± 0.2 ± 0.9 


17.1 ± 0.2 ± 1.3 


15.8 ± 0.4 ± 1.2 


1.0 


- 1.2 


14.5 ± 0.5 ± 1.2 


18.7 ± 0.5 ± 1.4 


15.1 ± 0.4 ± 1.2 


1.2 


- 1.4 


16.2 ± 0.4 ± 1.2 


18.9 ± 0.5 ± 1.5 


13.6 ± 1.1 ± 1.0 


1.4 


- 1.6 


17.8 ± 0.6 ± 1.3 


19.1 ± 0.7 ± 1.5 


12.6 ± 1.2 ± 0.9 


Long-track selection 


0.0 


-0.2 


0.8 ± 0.0 ± 0.2 


2.0 ± 0.1 ± 0.4 




0.2 


-0.4 


0.7 ± 0.1 ± 0.1 


2.0 ± 0.1 ± 0.4 




0.4 


-0.6 


1.2 ± 0.0 ± 0.2 


3.7 ± 0.1 ± 0.6 


1.3 ± 0.3 ± 0.2 


0.6 


-0.8 


1.9 ± 0.1 ± 0.3 


4.9 ± 0.1 ± 0.6 


2.9 ± 0.1 ± 0.4 


0.8 


- 1.0 


2.6 ± 0.1 ± 0.3 


5.6 ± 0.1 ± 0.7 


4.1 ± 0.5 ± 0.5 


1.0 


- 1.2 


2.8 ± 0.1 ± 0.3 


6.1 ± 0.5 ± 0.6 


4.3 ± 0.3 ± 0.4 


1.2 


- 1.4 


2.7 ± 0.2 ± 0.2 


5.7 ± 0.5 ± 0.5 


5.1 ± 0.6 ± 0.6 


1.4 


- 1.6 


2.8 ± 0.3 ± 0.2 


5.7 ± 0.5 ± 0.5 


5.4 ± 0.5 ± 0.5 



PV is required in the former case, so the systematic uncertainty on the PV reconstruction 
efficiency needs to be assessed. The simulation is found to be in good agreement with 
the data, but the analysis is more sensitive to the contribution from diffractive events. 
Secondly, the background level in the long-track selection is significantly lower than in 
the downstream-track selection, due to the PV requirement and the precise VELO mea- 
surements. Therefore it is possible to remove the minimum requirement on the K$ 
daughters in the long-track selection. This allows the extension of the analysis to two 
low j>t bins in the range 2.5 < y < 3.0, which are inaccessible to the downstream-track 
selection. The dominant systematic error for these two bins is from the large uncertainty 
on the tracking efficiency for the very low pr Kg daughters. 

The estimates of the total efficiencies e* rig,/scl x e^ cl are given in Table HI The various 
contributions to the uncertainties have been classified according to their correlations across 
bins, as shown in Table and added in quadrature. 
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Table 5: Sources of uncertainty on the K§ production cross-sections of Eq. ([I]), with relative 
values quoted for the downstream-track selection. A range of values means that the uncertainty 
was evaluated per bin of (px j y) phase space (with extreme values quoted) , while a single value 
indicates a global uncertainty assumed to be bin- independent. The different contributions are 
classified as uncorrelated or (at least partially) correlated across the different bins. 



Source of uncertainty 


uncorrelated 


correlated 


Yields N° bs 






- Data statistics 


5 - 25 % 




- Signal extraction 


1-5% 




- Beam-gas subtraction 




< 1% 


Efficiency correction (ef lg ^ sel e| el ) _1 






- MC statistics 


1-5% 




- Track finding 




6-17% 


- Selection 




4% 


- Trigger 




2% 


- pt and y shape within bin 


- 20 % 




- Diffraction modelling 




0-1% 


- Non-prompt contamination 




< 1% 


- Material interactions 




< 1% 


Normalization (L int ) _1 






- Bunch currents 




12% 


- Beam widths 




5% 


- Beam positions 




3% 


- Beam angles 




1% 


Sum in quadrature 


6 - 28 % 


16 - 23 % 



7 Results and discussion 

The cross-sections defined in Eq. ([1]) are evaluated separately for both the downstream- 
and long-track selections. In every phase-space bin, the two sets of results are found to 
be consistent with each other. The relative uncertainties on the measurement for the 
downstream-track selection are summarized in Table Since the downstream- and long- 
track results are not statistically independent, and since the downstream-track selection 
contains already most of the statistical power in bins where a measurement is possible, 
the measurements are not combined. The final results, listed in Table [6j are taken from 
the downstream-track selection, except in the two lowest px bins for 2.5 < y < 3.0 where 
they are taken from the long-track selection. 

The corresponding differential cross-sections are shown in Fig. as function of trans- 
verse momentum for the three different rapidity bins. They include both non-diffractive 
and diffractive prompt K§ production, and are compared with three different sets of 
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Table 6: Prompt K§ production cross-section (in fib) measured in bins of transverse momentum 
Pt and rapidity y, as denned in Eq. ([1]). The first quoted error is the statistical uncertainty, the 
second error is the uncorrelated systematic uncertainty, and the third error is the systematic 
uncertainty correlated across bins. 



Pt [GeV/c] 




2.5 < 


y <3.0 






3.0 < 


y <3.5 






3.5 < 


y < 4.0 




0.0-0.2 


294 


± 80 


± 38 


± 


90 


316 


±43 


±44 


± 


72 


196 


±39 


±39 


± 


38 


0.2-0.4 


649 


±133 


±136 


±1 


83 


562 


±42 


±22 


±101 


571 


±42 


±25 


±114 


0.4-0.6 


618 


± 63 


± 66 


± 


97 


534 


±30 


±12 


± 


86 


477 


±26 


±14 


± 


77 


0.6-0.8 


401 


± 32 


± 18 


± 


64 


371 


±21 


± 9 


± 


59 


323 


±20 


± 9 


± 


51 


0.8 - 1.0 


232 


± 20 


± 4 


± 


37 


183 


±14 


± 6 


± 


29 


201 


±16 


± 6 


± 


33 


1.0 - 1.2 


115 


± 11 


± 4 


± 


18 


136 


±11 


± 3 


± 


22 


108 


±12 


± 5 


± 


17 


1.2 - 1.4 


85 


± 10 


± 3 


± 


14 


70 


± 8 


± 2 


± 


11 


35 


± 9 


± 3 


± 


6 


1.4 - 1.6 


46 


± 7 


± 2 


± 


7 


49 


± 6 


± 3 


± 


8 


23 


± 6 


± 2 


± 


4 



predictions, all obtained with the PYTHIA 6.4 generator 13j. These predictions are 
represented as histograms in Fig. \5\ and correspond to: 

• the LHCb settings^), which include only soft diffraction as described by PYTHIA 6.4 
(red solid histogram); 

• the LHCb settings where diffractive processes have been switched off (blue dotted 
histogram) ; 



the "Perugia 0" settings 17|, which exclude diffraction (green dashed histogram). 



The predictions agree reasonably well with the data, although they tend to underestimate 
(overestimate) the measured production in the highest (lowest) pr bins. 

Previous measurements of the prompt K§ cross-section in high-energy pp collisions, 
performed by UA5 0, UA1 and CDF 0f at different centre-of-mass energies and in 
different rapidity or pseudo-rapidity regions, have been published in the form of invariant 
differential cross-sections E d 3 a/d 3 p as a function of pr- We convert these into measure- 
ments of d 2 a/ (dprch/) by multiplication with 2irpT, and compare them with our results 
in Fig. El limiting the px range of previous measurements to 1.6 GeV/c. In this figure, 
LHCb results are shown for the rapidity range 2.5 < y < 4.0, obtained by averaging the 
results for the three separate y bins, assuming conservatively that the correlated system- 
atic uncertainties are 100% correlated. In general the agreement is reasonable, given the 
spread of centre-of-mass energies and the fact that the results are averaged over different 



3 We use PYTHIA 6.421, and include process types 11-13, 28, 53, 68, 91-95, 421-439, 461-479 
with non-default parameter values ckin(41)=3.0, mstp(2)=2, mstp(33)=3, mstp(128)=2, mstp(81)=21, 
mstp(82)=3, mstp(52)=2, mstp(51) =10042, parp(67) =1.0, parp(82) =4.28, parp(89) =14000, parp(90) 
=0.238, parp(85)=0.33, parp(86)=0.66, parp(91)=1.0, parp(149)=0.02, parp(150)=0.085, parj(ll)=0.5, 
parj(12)=0.4, parj(13)=0.79, parj(14)=0.0, parj(15)=0.018, parj(16)=0.054, p_arj(17)=0.131, mstj(26)=0 
parj(33)=0.4. The particle decay probabilities are computed using EvtGe: 
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Figure 5: Double-differential prompt K§ production cross-section in pp collisions at -y/s = 
0.9 TeV as a function of transverse momentum and rapidity y. The points represent LHCb 
data, with total uncertainties shown as vertical error bars and statistical uncertainties as tick 
marks on the bars. The histograms are predictions from different settings of the PYTHIA gen- 
erator (see text). The lower plots show the MC/data ratios, with the shaded band representing 
the uncertainty for one of these ratios, dominated by the uncertainty on the measurements (the 
relative uncertainties for the other ratios are similar). 



ranges in rapidity or pseudo-rapidity. The ability of LHCb to contribute measurements 
that extend the kinematic range towards high rapidities and very low is apparent. 



8 Conclusions 

Studies of prompt K§ production at y/s = 0.9 TeV have been presented, made with the 
LHCb detector using the first pp collisions delivered by the LHC during 2009. The data 
sample used corresponds to an integrated luminosity of6.8±1.0/xb -1 , a value which has 
been determined using measurements of the beam profiles that exploit the high precision 
of the VELO. This is the most precise determination of the luminosity for the 2009 LHC 
pilot run, only limited by the uncertainties on the beam intensity. 

The differential cross-section has been measured as a function of pi and y, over a range 
extending down to pr less than 0.2 GeV/c, and in the rapidity interval 2.5 < y < 4.0, a 
region that has not been explored in previous experiments at this energy. These results 
show reasonable consistency with expectations based on the PYTHIA 6.4 generator, and 
should provide valuable input for the future tuning of Monte Carlo generators. 
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Figure 6: Absolute measurements of the prompt i£g production cross-section as a function of 
transverse momentum pr> performed by the UA1 [1], UA5 0, CDF 0] and LHCb experiments, 
at different high-energy hadron colliders and in different rapidity (y) or pseudo-rapidity (jj) 



ranges. 
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